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A GEOPHYSICAL ATLAS 

FOR INTERPRETATION OF SATELLITE-DERIVED DATA 

Paul D. Lowman, Jr. 

Herbert V. Frey 


INTRODUCTION 

Our knowledge of the earth’s structure has increased tremendously in the last two decades, 
much more so than would be expected in the normal course of scientific progress. There are several 
reasons for this increase, one of them being the achievement of space flight, and specific all y the 
launching of earth-orbiting artificial satellites. These satellites, either directly or indirectly, have 
permitted investigations of the earth’s gravity field, crustal magnetism, plate movements, and crustal 
structure. We present here the first unified compilation of satellite-derived and conventional geo- 
physical data applicable to such investigations. 

This Atlas consists of a series of maps intended for qualitative intercomparison. Several of 
these, such as the maps showing global seismicity, are ground-based data necessary for interpretation 
of the satellite data. The maps have been drawn on a common projection (Van der Grinten) and 
with a common scale, except for those showing the polar regions, which use an orthographic pro- 
jection. Each map represents the result of many years* work by many people and immense compi- 
lations of data (much of it digital), which would be impractical to reference completely. However, 
the main publications used are listed in the bibliography. In addition, the reader may find useful the 
Directory of U.S. Data Repositories Supporting the International Geodynamics Project” (Report 
SB-14, World Data Center A, National Oceanic and Atmospheric Administration, Boulder, CO 
80803, 1978). 

The Atlas is intended as a tool for scientists and educators rather than as a research publication 
in itself. For this reason, the authors have included only brief accompanying texts for the maps, de- 
signed to aid the reader in their use. A brief summary of the content and organization of the Atlas 
follows. 

The first map (Plate 1 ), reproduced with the kind permission of the National Geographic Soci- 
ety, shows the world’s physiography, both continental and oceanic. This map was chosen as a topo- 
graphic base because of its detailed representation of physiographic features such as ridges, trenches, 
and transform faults, the fact that it does not split any large land masses, and its Van der Grinten 
projection. The Van der Grinten projection, discussed in detail in Section 1 , though neither equal- 
area nor conformal, permits the inclusion of very high latitudes with moderate distortion compared 
to the commonly-used Mercator projection. A Van der Grinten outline map with latitude and lon- 
gitude is also presented (Plate 2). The first geophysical maps show various aspects of the earth’s 
gravity field, a subject of fundamental importance and, historically, the first branch of geophysics to 
benefit from satellite-derived data. Next follow maps of the earth’s crustal magnetism as measured 
by the Polar Orbiting Geophysical Observatory satellites. These are followed by a series of seismic 
activity maps, which are of interest both by themselves and because they served as a major input to 
other maps in the Atlas. Volcanic activity of the past one million years is shown both on a special- 
ized map and later on a combined map of tectonic and volcanic activity. We then present maps of 
the earth’s rifts and sutures, representing in plate tectonic theory the beginning and end of ocean 
basin evolution. These are followed by three maps showing global tectonic and volcanic activity for 
the past one million years. Finally , a series of ten maps shows seismic wave velocity anomalies for a 
set of five depth regions, from the surface to the core-mantle boundary at 2900 km depth. 
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Since this is the first such compilation of satellite-derived and related conventional geophysical 
data, much improvement in content and format is doubtless possible; we stress that this is a prelim- 
inary version and that improved editions will follow. We believe, however, that even this first edition 
of the Atlas will be useful as a research tool, an educational resource, and as a demonstration of the 
great scientific value of artificial satellites. 
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SECTION 1 : BASE MAPS AND PROJECTIONS 
John O'Keefe, Allen Greenberg, and Gilbert Mead 


Two projections have been used for all of the maps in this Atlas. For reasons given in the In- 
troduction, most of the maps were plotted, all to the same scale, on the Van der Grinten projection 
CVan der Grinten, 1905). The remaining data were plotted using an orthographic polar projection. 
The mathematical characteristics of these two projections are given here for the benefit of those 
who may wish to plot additional data to the same scale and projection. 


VAN DER GRINTEN PROJECTION 

A latitude-longitude grid of the Van der Grinten projection is shown in Plate 2. This plate, as 
weUasmost of the plates in this Atlas using the Van der Grinten projection, extends only to latitude 
±80 . The outline of the entire earth, including the poles, is a perfect circle. The projection has true 
scale along the equator. It is neither conformal nor equal-area. The characteristics of this protection 
were reviewed by O’Keefe and Greenberg ( 1 977), who also derived analytic equations for the rectan- 
gular coordinates (x,y) of a point P whose longitude measured eastward is X and whose latitude is 6 
TTie final results are given here ; for the specific derivation and the description of the projection as 
given initially by Van der Grinten, see O'Keefe and Greenberg (1977). 

Suppose that the longitude of the central meridian of the map is X<, ; we then define 


9 "o 

* " 180 

-1 < 2 < 1 

h - JL 
b " 180 

- Vi < b < 


The longitude scale is linear along the equator. Meridians are arcs of circles passing through the 
appropriate point on the equator and both poles. Their radius r is given by 

r-p(t + j) lrl>2p 

where the scateof the map is defined by the value assigned to p, which is the map distance corre- 
jp^ding to 90 of longitude along the equator. The center of the circle is along the equator at the 



The last two equations can be combined to give 

r 2 * x £ + 4 p 

Thus the ± 1 80° meridians (2 * ± 1 ) define a circle of radius r * 2p. The 90° east meridian (2 = 0.5) 
is the arc of a circle of radius r x 2.5p centered at x M * -l,5p and crossing the equator at x * p. 


PRECEDING PAGE BLANK NOT. FfcMED 
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The parallels are circles whose radius, s, is given by 


where z is defined by: 


_ 2pcosz 
sin 2 a 


and 


z * — lb| “ 90 ° < z < 90 ° 

sin z - tan p 


tana 


cos z 


sin 2 p = 2 b -45° < p < 45° 

The parallels intersect the prime meridian at 

y = 2 p tan p 

Thus the center of the circle defining the parallel is at 

yj = s + 2 p tan p 

- The gwmetry for the northeast quadrant is given in the attached figure. The x - y coordinates 
of the point P can be obtained by solving triangle MJP: y coordinates 

MP = r 


PJ = s 



Xy 

tan ij 0 * - - 0 < 77 < 2 w 

yj 

The angle M is defined by 

s 2 » r 2 + A 2 -2rAcosM 

Then 

x * r sin +x M 
y * rcostf. 
where 

n * V 0 +M 


Y 
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The geometry and formulas for the other three quadrants are given in the O'Keefe and Greenbeig 
paper. 


ORTHOGRAPHIC PROJECTION 

The orthographic polar projection used in Plates 12, 13, 18, and 19 shows the earth’s hemi- 
sphere as it would appear if viewed from infinity. In polar coordinates: 

r = p cos 0 
0 =X 

for the North polar cap, where 0 and X are latitude and longitude and the scale factor p, true only at 
the center of the projection, is the map distance corresponding to 1 80/w degrees of latitude. The 
sign of 6 is reversed for the south pole. 


REFERENCES 

Van der Grinten, A. J., New circular projection of the whole earth's surface, Am. Jour. Set, 19, 
357-366, 1905. 

O’Keefe, J . A., and A. Greenberg, A note on the Van der Grinten projection of the whole earth onto 
a circular disk, A mer. Cartographer, 4, 127-132, 1977. 
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VAN DER GRINTEN GRID 


SECTION 2: SATELLITE-DERIVED GRAVITY MAPS 


James G. Marsh 


INTRODUCTION 

Mapping the earth’s gravity field and determining its shape (the geoid) were among the first 
major scientific applications of data from artificial satellites. As early as 1958, a new value for the 
flattening of the earth was derived from the tracking of Sputnik II. Since then, progressively more 
detailed models of the earth’s gravity field anu the geoid (or equipotential surface) have been con- 
structed, based on hundreds of thousands of measurements of satellite orbits. We present here maps 
constructed from the latest of these, the Goddard Earth Model 10 (GEM 10). 

A helpful review of the use of satellite data in mapping the earth’s gravity field has been pub- 
lished by King-Hele (1976). The general subject of the earth’s shape and gravity field is covered by, 
for example, Garland (1965). 


DATA SOURCES AND COMPILATION METHODS 

The GEM 10 (Lerch, et al., 1977) is the latest in a series of models developed by Goddard 
Space Flight Center. Even-numbered models (“combination solutions”) are based on satellite and 
surface gravity data, odd-numbered models on satellite data alone. The GEM 10 model consists of 
a set of spherical harmonic coefficients complete to degree and order 22, with selected higher degree 
terms out to degree 30. It was derived from a combination of over 840,000 observations of 30 sat- 
ellites, of which 200,000 are precision laser ranges on 9 satellites, and a global set of 1654 5°x 5° 
mean free-air surface gravity anomalies. 

Three gravity anomaly maps are presented : one showing anomalies corresponding to the com- 
plete model, and two showing anomalies corresponding to degree and order 0-12 and 13-22. These 
maps are free-air anomaly maps, as contrasted with Bouguer maps. Among the major features shown 
by the full field map are: a major low south of India; a broad negative anomaly over the northwest 
Pacific; negative anomalies over Hudson Bay and Antarctica; positive anomalies over the Andes and 
Himalayas. The (0-12) map shows the long wavelength portion of the field. The India low is very 
prominent, and a positive ring clearly depicted around the Pacific basin closely matches plate bound- 
aries. The (13-22) map is primarily characterized by features with half wavelengths of about 1000 
km; some larger features shown on the (0-12) map have accordingly been removed. It has been pos- 
tulated (Marsh and Marsh, 1976) that the anomalies shown on the (13-22) map in the eastern Pacific 
may be indicative of convective rolls beneath the Pacific plates. The Mid-Atlantic ridge and the hot 
spot locations show no characteristic anomaly patterns. However, a prominent positive anomaly is 
apparent along the east African rift system. 

The fourth map presented here shows a detailed gravimetric geoid computed by combining the 
GEM 10 potential field model and a set of 38,406 1° x 1° mean free-air surface gravity anomalies. 
The long wavelength ('*' I500km> undulations were calculated using the GEM 10 gravity model and 
the short wavelength components were computed by applying Stokes' formula to 1° x 1° mean re- 
sidual anomalies. (The residual anomaly is the surface anomaly minus the satellite-derived anomaly.) 
Large geoidal lows (corresponding to negative gravity anomalies) are visible south of India, over Hud- 
son Bay, and in the Pacific Ocean west of Baja California. Significant highs (corresponding to posi- 
tive gravity anomalies) are visible over the Andes, south of Iceland, and over New Guinea. 
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SECTION 3: SATELLITE-DERIVED MAGNETIC ANOMALY MAPS 
H. V. Frey, R. A. Langel, and W. M. Davis 


INTRODUCTION 

Although not intended for such applications, the magnetic field measurements made by the 
Polar Orbiting Geophysical Observatory satellites (POGO-2, 4, and 6) were used by Regan et al. 
(1975) to construct a map of global crustal magnetism in the mid-latitudes. This map, the first such 
ever produced, revealed many anomalies, several of which have been confirmed by ground or aerial 
surveys. We present here modified versions of this map. 


DATA SOURCES AND COMPILATION METHODS 

This presentation represents a slight modification and expansion of the data presented by Regan 
et al. (1975), and is shown as contoured residual values after removal of a 1 3th degree and order 
main field (or core field). Two data sets were used, one with geomagnetic latitude limits of ±60°. 
the other with geographic latitude limits of ±50°. All the data were derived from POGO observations 
during magnetically quiet periods, excluding measurements during local daylight to avoid external 
field changes. These have been further treated in two ways to eliminate the pronounced north-south 
orbital striping of the first map published by Regan et al. These treatments are outlined briefly as 
follows; further details can be found in Langel (1974). 

Because typical signals measured at satellite altitudes (400-1 500 km for the POGO series) are 
comparable in strength to those due to time variations in external field sources (e,g.. ionospheric 
currents), these external sources must be removed from the data to determine the anomaly field due 
to internal sources. Contributions from distant external sources (above the ionosphere) and from 
long wavelength ionospheric sources can be modeled to first order by the low order terms in a spher- 
ical harmonic expansion (Langel and Sweeney, 1971 ; Cain and Davis, 1973). At middle geomagnetic 
latitudes (<60 ) the effects of fields introduced by ionospheric currents can be minimized by exclud- 
ing data from 0900 to 1550 local time. Restriction to <60° geomagnetic latitude also reduces con- 
tributions from polar and auroral zones. If these corrections successfully remove the effects of ex- 
ternal fields, the remaining signals are due to the Earth’s main field (i.e., the core field) and whatever 
lithospheric sources exist. The main field was modeled by Regan et al. by a 1 3th degree and order 
spherical harmonic expansion. Residuals from this field model are presumed to result from magnetic 
sources in the upper lithosphere. 

Despite the corrections described, the original POGO map showed strong north-south striping 
in the residuals, clearly related to the satellite orbital paths. Two approaches have been used to 
eliminate these distracting and spurious features, both based on 2° block averages rather than the 1° 
averages used by Regan et al. This coarser average represents very well the observed anomalies with 
little loss of spatial resolution. 

The first and simplest approach adopted was to further average out the orbital striping effects 
(which introduce high frequency noise along east-west lines) after removal of the external source 
model described above. Each value within a 2° block was adjusted by taking a weighted average with 
its nearest eight neighboring points; i.e., the data were smoothed in two dimensions using a nine- 
point weighted average. This clearly dampens the high-frequency noise introduced by the incom- 
plete removal of the long wavelength magnetospheric currents, and results in what will be referred to 
as a “smoothed” anomaly map. 
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The second approach was bused on Mayhew ’s (1977) study of individual POGO passes over 
the United States. Using a subset of the data with similar orbital altitudes, Mayhew found that 
individual profiles over the same region were not identical after correction for external effects. A 
long-wavelength residual remained in what should have been identical profiles; this was partly re- 
sponsible for the north-south striping in the original map. Mayhew found that an internally consis- 
tent data set could be produced by fitting linear or quadratic functions through the individual resid- 
ual profiles, and used that data set to generate an equivalent source representation of the anomaly 
field over the United States. This approach has been adopted here and applied to individual passes 
before the 2° averaging; the resulting plot will be referred to as a “filtered” anomaly map. 

Both approaches have shortcomings, and neither explicitly identifies the nature of the noise 
responsible for the striping. There are differences between the maps generated by the two tech- 
niques, which are particularly obvious for long wavelength features. In general, the “smoothing” 
procedure has probably failed to completely remove spurious external fields (e.g., magnetospheric 
effects), while the Mayhew filtering procedure may have removed these but in addition some of the 
lithospheric field features. 

Three anomaly maps are presented. The first (Plate 7), a smoothed anomaly map, represents 
2° block averages contoured at 2 gamma intervals between 60° geomagnetic latitude limits. This 
data set is sliehtly expanded over that used by Regan et al. (1975). The second map (Plate 8), a 
smoothed version, is similar to the first, but the anomalies are patterned with vertical and horizontal 
stripes reflecting the intensity of the anomalies. The zero contour has been eliminated for clarity. 
As in the first map, the contour interval is 2 gammas, but anomalies with absolute intensity greater 
than 6 gammas are shown only by 6 gamma stripes. The density of stripes in the pattern shows 
intensity of the anomalies, with horizontal stripes representing negative anomalies and vertical ones 
positive anomalies. This format readily displays the anomalies without the use of a color code such 
as that used by Regan et al. 

The third map (Plate 9) shows a filtered version of the anomaly data (patterned). This 
map is limited to <50° geographic latitude to accommodate restrictions in the Mayhew filtering 
technique. It is clear that most anomalies display only slight differences. However, there are some 
obvious changes, notably the anomalies in the Pacific west of Peru and Chile, in South Africa, and 
near Madagascar. Similar discrepancies for a small fraction of the anomalies occur if the quadratic 
filter recommended by Mayhew ( 1 977) is used instead of the linear trend removed here. 
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SCALAR MAGNETIC ANOMALIES 

Based on OGO 2. 4, 6 Data 
Filtered Version (Patterned) 

Contour Interval 21 
Van der Crinten Projection 
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SECTION 4: GLOBAL SEISMICITY MAPS 


M. K. Hutchinson 
P. D. Lowman 



compiled 

between «i U . wnuncmu ouiunes, me otner without. In addition, we have coi 

piled maps of the north and south polar regions, down to 40° latitude, to provide global coverage, 
fcach map can be compared directly with a corresponding tectonic activity map in this atlas. 


DATA SOURCES AND COMPILATION METHODS 

All the seismicity maps are computer plots derived from a data tape supplied by the National 
Geophysical and Solar-Terrestrial Data Center (NOAA) in Boulder, CO. This tape includes informa- 
tion such as data source, time of origin, location to three decimal places, focal depth, and associated 
phenomena if relevant. Seismic and other data 'apes available from NOAA are described in Report 
SE-14, published by World Data Center A for Solid Earth Geophysics, National Oceanic and Atmo- 
spheric Administration, Boulder, CO 80303. 

. P l ri0 £, waS chosen for P lottin 8 primarily because events from this period were re- 

corded by the World Wide Network of Standardized Seismographic Stations (WWNSSS) established 
in the early 1960 s. These data are more accurate than those from earlier times. Computer limita- 
tions restricted the number of data points that could be plotted, and the magnitude interval 4.5 to 
5.5 (generally M L ) was therefore chosen. As the maps demonstrate, despite the narrowness of this 
interval, earthquakes of this magnitude are quite numerous enough to outline zones and specific 
features of tectonic activity. One reason for this is that events larger than magnitude 5.5 have many 
aftershocks of lower magnitude. No discrimination was made on the basis of depth, with all earth- 
quakes down to 700 kilometers being plotted. Intermediate and deep focus events are abundant in 
subduction zones such as the Peru-Chile trench; outside such zones, almost all earthquakes are 
shallow-focus (0-70 km depth); (Tarr, 1 974). 

The map without continental outlines (Plate 10) includes 27,000 epicenters. Because of com- 
puter limitations, the inclusion of continental outlines required reduction in the number of points 
plotted; Plate 11 thus includes only 16,000 epicenters. To insure accurate delineation of seismic 
features, the reduction was accomplished by plotting the first 10 events in each 2 ° block, so that the 

net cartographic effect was simply a decrease in point density in the most active and best-defined 
seismic areas. 

These maps give a reasonably true and complete picture of the distribution of present seismic 
activity over the whole earth. Areas of little or no apparent seismicity in particular are probably 
areas in which tectonic activity at this time actually is minor. Antarctica is especially conspicuous 
in this regard; the near-absence of earthquakes from this continent is definitely not an effect of in- 
strument placement, seismometers having been operating there for some two decades. 

.. . I* 1 ? distribution of seismic activity along the oceanic ridges also deserves some comment First 
it is obvious that the activity varies drastically bothalong individual ridges and between ridges. What 
is not obvious is that the fastest spreading centers, in particular the East Pacific Rise, are not the 
most active ridges; seismicity apparently is greater along the Mid-Atlantic Ridge, whose spreading 
rate is much lower than that of the East Pacific Rise. Whether this is an artifact resulting from the 
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preceding page blank not filmed 



* he instrumental record - the narrow time interval plotted, or a real phenomenon, 


nQ7 ^f of the s ! 10I t time period represented by these maps, discussed at length by Allen 
( 76), is that some major faults known to be active are unrecognizable on the seismidtv Dlots For 
example, the South Atlas fault bounding the Atlas Mountains, which was responsible for the earth- 
quake that nearly destroyed the city of Agadir in 1961, exhibits only one event in the 1965-75 

Hi?* m , 3pS , by !!} emselves * hus d0 not for many areas give an accurate index of seismic risk 
?Pbt f , 7 ?S aCtm y ’ 7 ? e m u aps °J tectonic and volcanic activity presented elsewhere in this atlas 
(Plates 17-19) were therefore based on many criteria in addition to the seismicity maps. 
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SECTION 5: RECENT VOLCANIC ACTIVITY MAP 

Aron P. Trombka 
Paul D. Lowman 


This computer-plotted map (Plate 14) displays the geographical distribution of recent volcanic 
centers. Both currently-active volcanic centers and those active within the last one million years are 
included. 


DATA SOURCES AND COMPILATION METHODS 

The major reference for active volcanoes was the Catalogue of Active Volcanoes of the World 
by the International Association of Volcanology. Older volcanic centers were identified by using 
geologic and tectonic maps, especially those published by various national surveys. Recent informa- 
tion used to supplement outdated maps included the articles “Makran of Iran and Pakistan as an 
Active Arc System” by G. Farhoudi and D. E. Karig and “Concealed Mesozoic-Cenozoic Alpine 
Himalayan Geosyncline and its Petroleum Possibilities” by S. K. Achanya and K. K. Ray. For some 
isolated structures Landsat images were used to locate the features. 

The symbol (+) on the map represents one of these features: 

- All volcanic activity (central or fissure) listed as currently active by the International Asso- 
ciation of Volcanology or by other sources. 

- Any activity recorded (and confirmed by consulting a geologic map) during historical time. 

- An uneroded structure listed as “volcano” (A physical structure, e.g., a volcanic cone, that 
remains intact is probably less than one million years old, and as such is included here if the 
geology of the region, as indicated on a geologic map, supports the presence of young vol- 
canic material). 

- Lava or tephra fields definitely dated as formed within the past one million years. 

- Submarine volcanoes found within tectonic plates (hot spots) or those that form volcanic 
islands; however, all submarine volcanic activity associated with the mid-oceanic ridges and 
seafloor spreading has been omitted from this map. 

Any individual (+) indicates only the geographical point corresponding to the volcanic center 
and not the eruptive history of that center, i.e., repetitive activity from the same feature is plotted 
as a single (+). 


MAP LIMITATIONS 

Although this map is a complete compilation of all known active volcanoes, information for 
other features is less complete. Volcanic features other than volcanoes (e.g., solfatara fields, 
geysers), even while currently active, are in many instances never documented. 

Inactive volcanic centers present additional difficulties. Exact dating of volcanic features be- 
comes important when using such a brief time frame of one million years. Most of the younger 
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inactive volcanic center have P— ^ V 

ters with ages slightly less than one ™ lh °" J.^„ J nteR definitely dated as active within one mil- 
three million years. This map includes o y caseg where the exact location of older volcanic 
lion years and so is probably mcomplet . SV mbol (+) is plotted at the probable location, 

centers within a voicanic field is <^"e the symW M» discrepancy being 

this serves as an adequate representau™ on a map ot mis 
approximately three-quarters of a degree. 

Submarine volcanism is itself a f “ ly ^^"’^o^CothS man thoaScitted with sea- 
of specific events. Therefore, within the last one million years, 

noor spreading at mid^cean.c ndg l tuve been ^ed ^ activity has bee „ omitted from 

For both submarine and continental cases, reponea u 

this map. 
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SECTION 6: TECTONIC BOUNDARIES: RIFT AND SUTURE MAPS 

Herbert V. Frey 


INTRODUCTION 


'Die two maps (Plates 1 5 and 16) presented here were constructed by the Albanv Glohal TW 

siry ot New York at Albany). They are the only ones mtended to show major tectonic elements as 

HpW 1 ' a “ ap tf ropr ! a ! e t0 “ lclude specialized tectonic maps to aid correlative study of the satellite- 
d nved geophysical data, although it is not obvious just what types of maps will be m< st helpful 
Rifts and sutures are major tectonic elements representing, in plate tectonic theoiT The ooS and 
closing, respectively, of ocean basins, and they may delineate boundaries of lithospheric properties 

dtfto ,hUS rePreSent “ » fmt ste » f ” tacorpoSrfpKS 


DATA SOURCES AND COMPILATION METHODS 

The criteria used to identify rifts and sutures are those adopted by The Albany Global Tec- 

In TrSiSn U ^ aS SpCC,fied * n seve u ral Publications (Burke, 1977; Burke et al., 1977; Dewey,1977) 
If th» w m ^Ps are described and documented in a contractor report, “Rifts and Sutures 

NASA Sf ac «NA S a 5. e 24S COddard C «"’ ^ ~ pre“™der 

!ithn^S”i! *** extensional ruptures penetrating a substantial thickness of the 

lrthosptiere. A i ch ’ they represent major boundaries in lithospheric properties. This is espe- 
200 m!ninn fOr th ® continental nfts, which are the only ones preserved that a£j older than aboSt 
200 million years (the oceanic nfts are eventually destroyed by plate movements). Perhaps the best 
nown intra-continental rift system is that in East Africa, which includes the type Gregory rift 

ZZESftr SyS !»T S inC l Ude the Lake Baikal rift ’ the 05510 md RWne grabens, an® the Rio 
Graide nft. These and lesser-known nfts are mapped and described in the above-named contractor 


tonic^ctivkv mL d rthT P T if s Preadmg centers, and except for those now active (see the tec- 
), i h ? ftS Sh °'T on tlus map must be those which failed or whose development 
:rT^i bef TTr eIOpin8 .“! toa contmental mar gin. Some features displayed here are aulaco- 
KJ./ 1 QT 7 ?“ Th d ^ Wlt , h th j ck sedimentary sequences striking at high angles into fold belts 

by *nitaeS cSSSJT ' !ubse< ' uenUy f,Ued wi,h « d ™ ents “ d ifled 

. fv!l U L al co , n f < i u ® nce ° f P la | e tectonic theory is that continental blocks and other tectonic 
houl d eeibde. Complete closure of an ocean through subduction of the oceanic litho- 
sphere stops when contmental margins come into contact. The buoyancy of continental crust ore- 
S? 1 ? destruction through subduction. The weld where two previously-separated continental 
blocks join is often called, in plate tectonic theory, a suture. Sutures are seldom lines, but more 
often wide and complex zones displaying a variety of high-strain structures (Dewey, 1977). Ideally, 


•The Albany Global Tectonics Group includes Kevin Burke, 
K. D. Nelson, A. M. C. Sengor, and J. Stroup. 


Louise Delano, J. F. Dewey, A. Edelstein, W. S. F. Kidd, 
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TECTONIC BOUNDARIES SUTURES 

At Interpret*! by Burke et. al. (1978) 


SECTION 7: GLOBAL TECTONIC AND VOLCANIC ACTIVITY MAPS 

Paul D. Lowman, Jr. 


These maps are intended to show the geographic distribution of major tectonic and volcanic 
activity at the present time and for approximately the past one million years. There are three over- 
lapping maps, for the equatorial, Arctic, and Antarctic regions, which can serve as base maps for the 
qualitative interpretation of satellite data shown on other sheets of the Atlas. 


COMPILATION METHODS AND DATA SOURCES 

The equatorial tectonic activity map was drawn as an overlay on a 1 : 72,700,000 scale copy of 
the National Geographic Society’s “The Physical World” (1975). The actual delineation of tectonic 
and voicanic features was generally a compromise among several different representations, including 
published maps and papers listed in the bibliography. Plate boundaries, where well-defined, were 
drawn primarily from the epicenter patterns as plotted on the maps of global seismicity. Features 
included on the maps are described below. 


FEATURES DISPLAYED 
Active Ridges 


This category includes centers of sea-floor spreading and related continental rift zones. Activity 
is inferred from seismicity, but there is insufficient information in some areas to determine if appar- 
ently inactive ridge segments have been active within the past million years. An example of this 
problem is the segment of the East Pacific Rise at about 40°S, 1 10°W. Transform faults offsetting 
oceanic ridges arc highly generalized. Large oceanic fracture zones prominent on most maps, such 
as those of the eastern Pacific, are not shown here because they are now understood to be the inac- 
tive traces of transform faults. 

Subduction or Overthrust Zones 

Subduction and overthrust zones are both essentially reverse faults (though frequently associ- 
ated with fold belts), and they have accordingly been grouped together. This grouping may reflect a 
genetic relationship in many areas, such as the Himalayas, which are thought to represent a continent- 
continent subduction zone quite analogous to and continuous with the Java trench. Subduction 
zones have been located on the map over trenches representing their intersection with the surface of 
the earth. A strike-slip component of movement near subduction zones is shown only where prom- 
inent, as in the fold belts of Burma and Pakistan. 

Major Active Faults 


This ia a broad category for faults exhibiting good evidence (seismic or physiographic) of pres- 
ent or recent movement. Physiographic criteria, such as fresh scarps, were used by several authors 
whose work has been drawn upon, to infer fault activity; as pointed out by Allen (1975). such criteria 

are more reliable than the historic record of seismicity. Whether physiographic evidence of faulting 
necessarily indicates activity within the past one million years (the criterion for inclusion on this 
map) is an open question which is beyond the scope of this discussion; the reader is referred to lit- 
erature cited for individual structures. 
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Rifts or Zones of Normal Faulting 


This category includes areas of well-documented extensional tectonism not necessarily related to 
oceanic spreading centers. These features have been highly generalized in many areas, and some in- 
dividual rifts have been shown with exaggerated widths for clarity. Fault-controlled valleys of the 
Basin and Range province are shown only schematically, since it was impossible to represent them 
all at this scale. 

Young Volcanoes 


This category includes chiefly volcanic areas ot the central eruptive type, although a few 
(notably Iceland) have associated fissure eruptions. The problem of whether a volcano now inactive 
has erupted within the past one million years is for many areas a difficult one. It is obvious that this 
map and the “Recent Volcanic Activity” map (Plate 14) include many more volcanoes than those 
generally catalogued as “active,” though omitting many features shown as hot spots by Burke and 
Wilson (1976) (whose age cut-off was 10 million years). Inclusion in the one million year age cate- 
gory was largely on the basis of physiography. It appears, from study of the published literature 
(e.g., Hopkins, 1963) , that volcanoes in general survive as identifiable land forms for roughly a mil- 
lion years; i.e., if a feature is labeled “volcano” on a reliable map it is probably no older than this. 

It should be noted that active oceanic ridges are also the sites of igneous activity, although this 
is not shown explicitly on the map except for a few major volcanic centers such as Iceland. 

Unresolved Problems 

As a service to the user, some of the weaknesses of this map should be pointed out. A funda- 
mental inconsistency occurs between the seismicity and tectonic activity maps, in that a number of 
areas exhibiting seismicity are not represented on the tectonic activity map; examples include the 
Appalachians, the Urals, and the Canadian Arctic islands. This omission was caused by the difficulty 
of assigning the seismicity patterns to discrete tectonic features. In some areas, particularly the Ap- 
palachians, seismic activity may result from reactivation of Paleozoic structures. The problem will 
be addressed in future editions of the Atlas. 

A related problem is the treatment of “passive” continental margins, such as those of the Atlan- 
tic Ocean. The actual degree of activity over the last million years is quite unknown, and there is 
growing evidence of Tertiary or Quaternary deformation, for example, along the margin of eastern 
North America (Mixon and Newell, 1977) and of strong regional stress in this area (Sbar and Sykes, 
1973). This problem is a fundamental one whose solution must await more information about the 
nature of these features. 


DISCUSSION 

These tectonic activity maps fill a cartographic gap between conventional geologic maps and 
maps showing present seismic and volcanic activity. Although based partly on compilations of cur- 
rent activity, the use of geologic evidence of recent age permits extrapolation back about one million 
years, giving a much more comprehensive view of global tectonic and volcanic activity. An example 
of this is furnished by the many volcanic areas shown in the Alpine fold belt, very few of which are 
included in conventional compilations. Although most of these volcanoes are not now active, they 
are obviously related to the Alpine orogeny. 

The tectonic activity maps show most of the generally recognized tectonic plates and related 
features. However, they also show many structures generally omitted from conventional plate maps 
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such as that of Le Pichon et al. (1972), especially active structures located within plates or diffuse 
i P Jj? e Examples of these include the faults of western China, the Bastoand Range prov- 

ince, and the ‘trike-clip faults of Alaska. Omission of such features can obscure important relation- 

& « U965?' fundan,ent, ‘ l sym,nctr >' of "« Alpine fold belt in the Middle East discussed by 

In summary, the tectonic activity maps presented here provide a more realistic picture of crus- 
dynamic processes on a global basis than do conventional maps. They should be useful for geo- 
physical correlations as well as for tectonic research and education. g 


REFERENCES 

A|b any Global Tectonic Group, Rifts and Sutures of the World, NASA Contract Report NAS5-24094 
Goddard Space Flight Center, Greenbelt, Md., 238 p., 1978. ’ 

AHen i057 R i?7°5 ,08iCal Criteria ,or coating seismicity, Geol. Soc. America Bull., v. 86, p. 1041- 

AmC ’ Th f F, T r °, f th A e ° ceans (based studies by B. C. Heezen and 
10032 T974 40 ' 000 ’ 000 eqUat ° nal sca,e ’ Amencan Geographical Society, New York, N.Y. 

AndC n r; o R f "; SJhS’ n ngSeth / !? d J , G - l dater ’ 7116 ,nec hanisms of heat transfer through the 
floor of the Indian Ocean, J. Geophys. Res., v. 82, p. 3391-3409, 1977 

Atwater, T„ Implications of plate tectonics for the Cenozoic tectonic evolution of western North 
America, Geol. Soc. America Bull., v. 81, p. 3513-3536, 1970. 

Ben Avraham, Z.^and A. Nur, Slip rates and morphology of continental collision belts, Geology, 4, 

Barke i972 F ” A SprCad ‘ ng CCntrC in thc East Scotia Sca ’ Earth and planet - Sci. Lett., 1 5. 1 23-132, 

Bergli H. W and I. O. Norton, Prince Edward fracture zone and the evolution of the Mozambique 
Basin, J. Geophys. Res., v. 81, p. 5221, 1976. q 

BUrke i 97 6 C ’ and J ‘ T ‘ Wi,son ’ Hot spots on thc Earth ’s surface, Sci. American, 235. 46-57, August 

ChC Ti™ ’ l*975 ld W ' K ' Hamblin ’ P,, y si °g ra P hic Ma P of Earth, Burgess Pub. Co., Minneapolis, 

Churkin, M., Jr., Western boundary of the North American continental plate in Asia, Geo! Soc 
America Bull., v. 83, p. 1027-1036, 1972. 

Da,zi ^.;- D ' E A vo L uti ? n °[ the margins of the Scotia Sea, in The Geology of Continental Margins, 
edited by C. A. Burk and C. L. Drake, pp. 567-569, Springer-Verlag, New York, 1009 p., 1974. 

Dewey, J. F„ Plate tectonics, Sci. American. 226, 56-68, 1972. 


37 


Dewey, J. F., and J. M. Bird, Mountain belts and the new global tectonics. Jour. Geophys Res 
v. 75, p. 2625-2647, 1970. 

Dewey, J. F., W. C. Pitman, III, W. B. F. Ryan, and J. Bonnin, Plate tectonics and the evolution of 
the Alpine System, Geol. Soc. America Bull., v. 84, p. 3137-3180, 1973. 

Forsyth, D. W., Fault plane solutions and tectonics of the South Atlantic and Scotia Sea, J Geophvs 
Res., 80, 1429-1443, 1975. 

Gansser, A., Geology of the Himalayas, Interscience Publishers, New York, 289 p., 1964. 

Gutenberg, B., and C. F. Richter, Seismicity ot the Earth and Associated Phenomena, Princeton 
University Press, Princeton, NJ, 273 p„ 1949. 


Hilde, T. W. C., S. Uyeda, and L. Kroenke, Tectonic history of the western Pacific, in Geodynamics: 
Progress and Prospects, edited by C. L. Drake, pp. 1-15, Am. Geophys. Union, Washington, DC, 
238 p.,1976. 


Hinze, W. J., L. W. Braile, G. R Kellei, and E. G. Lidiak, A tectonic overview of the central midcon- 
tinent (preprint) 98 p„ 1977. 

Holmes, A., Principles of Physical Geology, Second Edition, Ronald Press, New York, 1288 p.. 1965. 

Hopkins, D. M., Geology of the Imruk Lake area, Seward Peninsula, Alaska, U.S. Geological Survey 
Bulletin 1 141-C, Contributions to General Geology in 1961, 1963. 

Isaacks, B., J. Oliver, and L. R. Sykes, Seismology and the new global tectonics, J Geophvs Res 
73, 5855-5899, 1968. 

Johnson, G. L.. and P. R. Vogt, Marine geology of Atlantic Ocean north of the Arctic Circle, in 
Arctic Geology, edited by M. G. Pitcher, Mem. 19, pp. 161-170. Am. Assoc. Petroleum Geo- 
logists, Tulsa, Oklahoma, 747 p„ 1973. 

Kaula, W. M., Global gravity and tectonics, in The Nature of the Solid Earth, edited by E. C. Rob- 
ertson, J. F. Hays, and L. Knopoff, pp. 385-405, McGraw-Hill, New York. 677 p., 1972. 

King, P. B., and G. J. Edmonston, Generalized Tectonic Map of North America, 1 : 1 5,000.000: Map 
1-688, U.S. Geological Survey, Reston, Va., 1972. 

Kirkham, R. M. Quaternary movements on the Golden fault, Colorado. Geology. 5. 689-692, 1977. 

Lathram, E. H., Nimbus IV view of the major structural features of Alaska. Science, / 75 1423-14'*7 
1972. 

Le Pichon. X. J. Erancheteau, and J. Bonnin, Plate Tectonics. Elsevier, New York. 300 p.. 1973. 

Eowman, P. D.. Jr., The Third Planet. Weltflugbild, Zurich, Switzerland. 185 p.. 1972. 

McKenzie, D. P„ Speculations on the consequences and causes of plate motions, Geophvs Jour. 

Royal Astronomical Soc., 18, 1-32, 1969. 

McKenzie, D. P., and W. J. Morgan, Evolution of triple junctions. Mature, 224. 1 25-1 33, 1969. 


38 


Meyerhoff, A. A.. Origin of Arctic and North Atlantic Oceans, in Arctic Geology, edited by M. G. 
Pitcher, Mem. 19, pp. 562-583, Am. Assoc. Petroleum Geologists, Tulsa, Oklahoma, 747 p., 


Mixon, R. B., and W. L. Newell, Staflord fault system: structures documenting Cretaceous and 
Tertiary deformation along the fall line in northeastern Virginia, Geology, 5, 437-440, 1977. 

Molnar, P„ and P. Tapponier, Relation of the tectonics of eastern China to the India-Eurasia colli- 
sion: application of slip-line field theory to large-scale continental tectonics. Geology 5 2P- 
216, 1977. ' 

Muehlberger. W. R., and A. W. Ritchie, Caribbcan-Americas plate boundary in Guatemala and 
southern Mexico as seen on Skylab IV orbital photography, Geology, 3, 232-235, 1975. 

O Keefe, J. A., and A. Greenberg, A Note on the Van der Grinten Projection of the Whole Earth 
Onto a Circular Disk, The American Cartographer, v. 4, no. 2, p. 1 27-1 32, 1977. 

Powell, C. McA., and P. J. Conaghan, Tectonic models of the Tibetan plateau. Geology, 3, 727-731, 


Prostka, H. J., and S. S. Oriel, Genetic models for Snake River plain, Idaho (abstract). Abstracts With 
Programs, 1975 Annual Meeting, Geol. Soc. America, v. 7, no. 7, p. 1236, 1975. 

Ray, K. K., and S. K. Acharyya, Concealed Mesozoic-Cenozoic Alpine-Himalayas geosyncline and 
its petroleum possibilities,^/??. Assoc. Petroleum Geologists Bull., 60, 794-808, 1976. 

Rutten, M. G., The Geology of Western Europe, Elsevier, New York, 520 p., 1969. 

Sbar. M. L., and L. R. Sykes, Contemporary compressive stress and seismicity in eastern North 
America: an example of intra-plate tectonics, Geol. Soc. Amer. Bull., 84 , 1861-1882, 1973. 

Sclater, J. G.,C Bowin, R. Hey, H. Hoskins, J. Peirce, J. Phillips, and C. Tapscott, The Bouvet triple 
junction,/ Geophys. Res., 81, 1857-1869, 1976. 

Simpson, R. W., and A. Cox, Paleomagnetic evidence for tectonic rotation of the Oregon Coast 
Range, Geology, 5, 585-589, 1977. 

Smith, P. S., Mineral resources of the Lake Clark-Iditarod region, in Mineral Resources of Alaska 
Reports on Progress of Investigations in ’914, edited by A. H. Brooks, pp. 247-271 , U S 
Geological Syrvey Bulletin 622, 1915 

Stein, S„ and E. A. Okal, Seismicity and tectonics of the Ninetyeast Ridge area: evidence for inter- 
nal deformation of the Indian plate, / Geophys. Res., 83, 2233-2245, 1978. 

Swiss Reinsurance Company, Atlas on Seismicity and Volcanism, Kummerly-Frey Pub. Co., Zurich, 


Tapponier, P., and P. Molnar, Active faulting and tectonics in China, / Geophys Res 82 "’905- 
2930. 1977. “ 

White, R. S., and K. Klitgord, Sediment deformation and plate tectonics in the Gulf of Oman 
Earth and Planet. Sci. Lett., 32. 199-209, 1976. 


39 


ORIGINAL PAGE 13 
Of POOR QUALITY 


ORIGINAL PAGE IS 
OF POOR QUALITY 







Pacific 




1^.0 N° h S 

y , ,) 


Plate qj 


O ^ ^ < .j 

'r f7A^ 


Wj? Pocific 
Jj Plate 


Equatorial Scale 


Cocos 5 
• Plate 


**■ H Naze a 


Plate ?/' 




t- *£• ^ 




Phil*)pr>e ^ 
’ Plate 


to, ^ 



Plate 


Antarctic 

Plate 


Amelncan » 

^Pla.e flf 


Sandwich 

Pli-:«, 


V f'V 

U rucan 1 

/ '<i J • y\ 

Plfate f y 


ScotiaPlo|* [J j 


I Amfrafcon 
i Plate 


GLOBALTECTONIC AND VOLCANIC ACTIVITY 
OF THE LAST ONE MILLION YEARS 

Goddard Space Flight Center 
1978 

Bated on ’The ^Nvcol 'Mxld." copy right 1975 by fhe Nanonoi 
Geograph* Socwty 
Von de. Grr**n Protector 

Paul D lowmon Jr 




Antarctic Plate 




Active ridges ond coni men toi extensions, minor faults 
f y _ generalized 

Subdue tion or overthrust zone, teeth on upper block 
Hajor octive foult or foult zone, dashed where nolure 

^ v octnntjr uncertain 

^ Rift or zone of normal faulting, teeth on downthrown block 

Young volcanos (octive wiffun last I million years, generalized) 








SECTION 8: SEISMIC VELOCITY ANOMALY MAPS 


Barbara E. Lowrey 


Long wavelength gravity anomalies such as those shown in Plate 4 reflect the mass distribution, 
or variation in density, in the deep interior of the Earth. To aid in understanding the causes of 
gravity anomalies, we have constructed a series of maps (Plates 20-29) showing anomalies in the ve- 
locities of seismic P-waves. 

The velocity of seismic waves has for many years been a primary source of information on the 
structure of the earth’s interior. Seismic wave travel times have been employed by many workers to 
determine the radial variation of density within the earth, assuming a spherically symmetric earth. 
However, it has for some time been recognized that there are lateral heterogeneities in the mantle. 

A major investigation has been undertaken by Dzicwonski, Hager, and O’Connell (1977) to deter- 
mine the feasibility of recovering the lateral velocity variations, i.e., the heterogeneities, in the man- 
tle on a global basis by computing the differences between observed P-wave travel times and the 
standard Jeffreys-Bullen travel time tables. Dzicwonski et al. constructed a spherical harmonic ex- 
pansion of the velocity anomalies, and published a series of Mercator projection maps showing the 
anomalies in successively deeper shells. 

Because the conversion of the recovered data into spherical harmonics truncated at degree and 
order 3 may lose numerical information, we have constructed contour maps directly from the re- 
covered data, using the Van der Grinten projection common to the maps of this Atlas. 


DATA SOURCES AND COMPILATION METHODS 

The work of Dziewonski et al. was based on 12,882 earthquakes for the years 1964-1970 pub- 
lished in the Bulletin of the International Seismological Centre. Compression wave (P and PKP) 
times were used, being the most accurate and numerous. Residuals were computed by differencing 
the observed travel times with the Jeffreys-Bullen travel time tables. The mantle was divided into 
five concentric shells, with depths of 0-670, 670-1 100, 1 100-1500, 1500-2200, and 2200-2886 kilo- 
meters, and 36° wide latitudinal and 60° wide longitudinal zones. Two cases were investigated, the 
second involving rotation of meridional boundaries 30° to the west to see if the results had been 
biased by the longitudinal boundaries selected. 

Dziewonski et al. reported that their results were unreliable for depths less than 1 100 km, prob- 
ably because lateral heterogeneities in this shell are generally smaller than 5000 km in lateral dimen- 
sions. However, they found a statistically significant correlation between long wavelength gravity 
anomalies observed at the surface and those computed from velocity anomalies deeper than 1 100 km, 
assuming a proportionality between velocity and density perturbations. Their interpretation of these 
results will not be repeated here. 

In the paper by Dziewonski et al., contour plots were constructed after several steps of numer- 
ical processing. First, spherical harmonic coefficients were computed according to two methods, 
numerical integration and least squares solution. Thus four sets of coefficients were computed for 
the two cases by two methods. The coefficients from the two cases were then combined for the 
two methods. Contours were constructed for each of the two methods and plotted together for 
comparison. 

The contours on our maps were constructed for the two cases directly from the original data 
(velocity anomalies) tabulated by Dziewonski, without additional numerical treatment. This has 
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several advantages: ( 1 ) loss of information by truncation of the spherical harmonics is avoided, 

, “"certainties due to lack of data or poor resolution are mapped directly, rather than as artifacts 
elsewhere onthe map, and (3) the tendency of spherical harmonics to produce non-existent anomalies 
is avoided. The method of contouring the original data is also useful in assessing the validity and 
degree of knowledge of a feature. Since the original study used large blocks (6 (f x 60°) to discretize 
the data, the geographic resolution of the data is about 30°. By comparing maps of a given depth 
region for the two cases, the size and strength of individual features can be assessed. In the upper 
two regions, there is little similarity between maps for the two cases, i.e., it is difficult to identify 
features common to both cases, probably because the resolution of the data is too coarse. However 
in the lower three regions there is a strong correlation between the two cases, although features may 
be offset as much as 30 because of the spatial resolution. 

°u ' treatme , nt ^/roduces a certain degree of complexity and ambiguity to interpretation of the 
data. However i it is clear that the general features of the maps by Dziewonski e* al. were reproduced 
at least for the deep regions of the mantle. Our maps are intended not to replace those of Dziewonski 
et al but to complement them by contouring the original values, thus providing a better idea of data 

onhe y iithmnf nbU h 10 m I" f dd r 10n ’ qualitative correlation between deep mantle structure and that 
In. lithosphere should be facilitated by the use of a scale and projection common to other geo- 
physical maps, particularly gravity, in this Atlas. oiner geo- 
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SEISMIC VELOCITY ANOMALIES (m/s) 
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Based on Data from Dziewonski. 
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Van der Grinten Projection 
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SEISMIC VELOCITY ANOMALIES (m/s) 
Case 1: Depths 670 to 1100 km 
Based on Data from Dziewonski. 

Hager, and O'Connell 
Van der Grinten Projection 

F T Heuring Goddard Space Right Center 
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Case Z Depths 1100 to 1500 km 
Based on Data from Dziewonski, 


Hager, and O'Connell 
Van der Grinten Projection 
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SEISMIC VELOCITY ANOMALIES (m/s) 
Case 1: Depths 1500 to 2200 km 
Based on Data from Dziewonski, 

Hager, and O'Connell 
Van der Grinten Projection 
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SEISMIC VELOCITY ANOMALIES (m/s) 
Case 1: Depths 2200 to 2886 km 
Based on Data from Dziewonski, 

Hager, and O'Connell 
Van der Grinten Projection 

r. T. Heuring Goddard Space Right Center 
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Case 2: Depths 0 to 670 km 
Based on Data from Dziewonski, 

Hager, and O'Connell 
Van der Grinten Projection 
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SEISMIC VELOCITY ANOMALIES (m/s) 
Case 2: Depths 1100 to 1500 km 
Based on Data from Dziewonski, 

Hager, and O'Conoell 
Van der Grinten Projection 

F. T. Heuring Goddard Space Right Center 
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